Magnetic subband structure of an electron on a triangular lattice 
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The effects of anisotropic nearest-neighbor (NN) and isotropic next-nearest-neighbor (NNN) 
hopping integrals on the magnetic subband structure of an electron on a triangular lattice are 
studied within the tight-binding approximation. A generalized Harper equation that includes both 
the NN and the NNN hopping terms is first derived. Then, the effects of anisotropic NN hopping 
integrals are examined and the results are compared in detail with those of an existing literature 
[Phys. Rev. B 56, 3787 (1997)]. It is found that the hopping anisotropy generically leads to the 
occurrence of band broadening and gap closing. The effects of isotropic NNN hopping integrals are 
next examined, and it is found that introducing the NNN hopping integral changes considerably 
the magnetic subband structure; band broadening, gap closing, gap reopening, and band crossing 
occur depending on the strength of the isotropic NNN hopping integral. Symmetries of the magnetic 
subband structures with and without both the hopping anisotropy and the isotropic NNN hopping 
integrals are also discussed. 

PACS numbers: 71.28.+d, 71.20.-b, 73.20.Dx, 71.45.Gm 
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I. INTRODUCTION 

The problem of a two-dimensional (2D) Bloch elec- 
tron under a perpendicular magnetic field has attracted 
much interest for several decades, since not- only it ex- 
hibits extremely rich energy band structuresuQ but also 
it is closely related with various phenomena such as the 
quantum Hall cffqetjj the flux-state model for high-T c 
superconductivity,!^! and the mean-field transition tem- 
perature of superconducting networks or Josephson junc- 
tion arrays. L3 Recently, by means of the remarkable ad- 
vent of nanofabrication techniques, studies of finding the 
indication of the peculiar band structure and its effect 
on the transport and .optical properties have been also 
extensively performedlLirtLi 

A traditional method to solve this problem is to start 
with energy dispersion e(k) without a magnetic field and 
then make the Peierls substitution k — > (p + eA)/h to 
construct an effective single-band Hamiltonian 



H = e[{p + eA)/h], 



(1) 



where p is the momentum operator and A is a vector 
potential. An equivalent way to this method is to start 
directly with the tight-binding Hamiltonian given by 



H 



E 



t; j e 



iO, 



(2) 



where £y is the hopping integral between the sites i and 
is a state of an atomiclike orbital centered at the 



J, 

site i, and Oij is the magnetic phase factor defined by 



2n f 3 a ,r 

— A-dL 



(3) 



So far, a number of works have focused on the energy 
spectrum of a lattice with a square symmetry, and the 
magnetic subband structure of the lattice with isotropic 
hopping integrals between nearest-neighbcffj/JMN) sites is 
well known as the Hofstadter's butterflyoca The effect 
of the anisotropy in the NN hopping integrals on the 
magnetic subband structure was also studied, and it was 
found that introducing the hopping anisotropy leads to. 
the occurrence of band broadening and gap closing. EjEij 
Furthermore, the effect of next-nearest-neighbor (NNN) 
hopping integrals on the magnetic subband structure was 
also studiedjij and introducing the NNN hopping inte- 
grals was found to play a key role in removing the de- 
generacy which appears at the band center of the energy 
spectrum under isotropic NN hopping integrals. 

Along with the case of the square lattice, magnetic sub- 
band structures, pf a-triangular lattice have been inten- 
sively studiedJi3'E2rE-il and it was found that the lattice 
with isotropic NN hopping integrals exhibits a recursive 
magnetic subband structure similarly to the case p£ the 
square lattice with isotropic NN hopping integrals .Ell Re- 
cently, the effect of the anisotropy in the NN hopping 
integrals on the magnetiC|-s.ubband structure was stud- 
ied by Gumbs and FeketeBJ In studying the effect, they 
applied the Peierls substitution for the energy dispersion 
given by 



e(k) = 2 |t a cos(fc;c a ) + h cos \^{k x + VSk y ) ' 



+t c cos 



(k x - V3ky)^ | 



(4) 



6q = he/ e being the magnetic flux quantum. 



to construct a Hamiltonian matrix. And, by diagonal- 
izing the resultant matrix, they obtained magnetic sub- 
band structures for a few sets of (t a , i&, t c ) and presented 
relevant arguments. However, unfortunately, most of the 
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results and arguments presented by them are erroneous 
ones that arise from a mistake in deriving the Hamilto- 
nian, as we shall show below. We re-examine the same 
problem as that investigated by Gumbs and Fekete to put 
the correct formalism and numerical results, which is the 
first end of this paper. The second end of this paper is 
to examine the effect of isotropic NNN hopping integrals 
on the magnetic subband structure. To our knowledge, 
few work has been devoted to this problem. However, 
noting that the triangular lattice with NNN hopping in- 
tegrals is topologically equivalent to the square lattice 
with NNN hopping integrals, it may be interesting to 
examine this problem and compare the results with the 
case of the square lattice. We will show that the mag- 
netic subband structure changes considerably depending 
on the strength of the isotropic NNN hopping integrals. 

This paper is organized as follows: In Sec. II, we de- 
rive a generalized Harper equation that includes NNN 
as well as NN hopping integrals. In Sec. Ill, we first 
point out the mistake and the misleading arguments pre- 
sented in Ref. 24, and then present correct numerical re- 
sults on the effects of both anisotropic NN and isotropic 
NNN hopping integrals on the magnetic subband struc- 
ture. We also discuss the symmetries of the magnetic 
subband structure obtained from the generalized Harper 
equation. Finally, Sec. IV is devoted to a brief summary. 



and (|5|), the tight-binding equation, H^ m n = E^ mn , 
can be written as 

+t 6 [e i »*('»-V2)^ m _ lin _ 1 + e -^ m+1 / 2 ^ m+hn+1 ] 

+tc [ e ^(m+l/2)^ m+ln _ 1 + e ^H m -l/2)^ m _ in+i] 

+.f> r p iT^(m-3/ 2 )fl/ „ , -^(m+S^j, I 

'^a[^ ^ m — 3,n — 1 i ( - ^m+3,n+lj 

+t' c {e l7rHm+3/2) *m+3,n-l + e-^ m - 3 / 2 >* m _ 3 ,„ + l] 

+t ' b ^m^ mn _ 2 + e - l2 ^ m * m ,„ +2 ), (6) 

where t a is the NN hopping integral along the x direction, 
t b ( c ) is the NN hopping integral along the direction which 
makes 7r/3 (27r/3) with the x direction, and t' a , b c ^ is the 
NNN hopping integral along the direction which makes 
ir/6 (n/2, 57r/6) with the x direction, respectively. 

Since Eq. (0) is invariant under the transformation 
y — > y + \/3a, ^> m ,n can be written as 

*m,„ = e^^Vw (7) 

Thus, using Eqs. @ and (0), we obtain a kind of gener- 
alized Harper equation as follows: 

Elp m = C* m _tf\) m -?, + t a -0m-2 + S^.^ro-l + A m ip m 

+B m ip m+ i + t a ip m+ 2 + C m+ i-0 m+ 3, (8) 



II. THE TIGHT-BIDING EQUATION 



where 



We consider an electron on a 2D triangular lattice with 
a hexagonal symmetry in the presence of a uniform per- 
pendicular magnetic field B = Bz. Under the Landau 
gauge, the vector potential is given by A = (0, Bx, 0) 
and 9ij is given by 



r o, j 

±7T0(m±l/2), j 

±7n^(m =F 1/2), j 

±tt0(to±3/2), j 

±7r^(mT3/2), j 

±27r0m, j 



[m ± 2, n) 
(m ± 1, n ± 1) 
(m =p 1, n ± 1) 
(m ± 3, n ± 1) 
(m =p 3, n ± 1) 
(m, n ± 2) 



(5) 



where i = (m,n). Here, (m, n) labels the lattice points, 
i.e., (x,y) = (mb,nc), where b = a/2, c = y/3a/2, and a 
is the lattice constant. And, 4> — 2^/3Bb 2 /4>o is the mag- 
netic flux through the unit cell. By means of Eqs. (H) 



A m = 2t' b cos(29 m - ncj)), 



B =t e ien 

v - y rr). 



+ t b e 
+ t'e 



-i0„ 



(9) 



with 6 m = TT(j)(m — 1/2) — k y c. Denoting cj) = p/q with 
relative primes p and q, it can be easily checked that 

A m , B m +M — B m , C m +M — C m , (10) 

where M = q (2q) for an even (odd) p. Thus, m in 
Eq. (||) satisfies the condition 1 < m < q (I < m < 2q) 
for an even (odd) p, and the Bloch condition along the x 
direction can be written as 



ipm+M 



e lk * Mb iP r , 



(11) 



Using Eqs. (|8j) and (|ll|), we obtain the eigenvalue equa- 
tion, A* = _E<3>, where *t = (-0 1; ip 2 , ■ • • , V'm) and 
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with S — k x qa/2 (k x qa) for an even (odd) p. 

6 rn — i — i — i — i — i — i — i — r 
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FIG. 1. Energy eigenvalues versus (f> for a triangular lattice 
with isotropic NN hopping integrals (t„ = tb = tc = 1). Cal- 
culations are performed for = p/199 with 1 < p < 198, and 
the central value of k in the MBZ is taken into account. 



III. NUMERICAL RESULTS AND DISCUSSION 

Before presenting our numerical results, we would like 
to point out the mistake made by Gumbs and Fekete 
in Ref. 24. In converting Eq. (Q) into Eq. (jl]) [i.e., 
Eq. (2) of Ref. 24], they made a mistake in operating 
e ±i( Px ±V3eBx)a/2h tQ ^ n _ According to them, the re- 
sults of the operations are given by 

±i(p :c +V3eBx)a/2h-fl l ±i7r0m,T. 

±i(p J: -V3eBx)a/2h^j _ „T^0™vI/ , fl ^ 

However, these are incorrect since the fundamental com- 
mutation relation = ih is ignored in the opera- 
tion. The correct results-taking into account this relation 
should be expressed aso 

±i(p x + V3eBx)a/2h x u _ ±1^^1/2)^ 
±i(px-V3eBx)a/2h. k T ) _ ^47r0(m=Fl/2) fl> (I 4\ 

We can easily check that the same equation as Eq. (|^) 
with t' a = t' b = t' c = can be obtained by the method of 
the Peierls substitution if Eq. (|l4|) instead of Eq. ( |l3| ) is 
applied. In fact, Eq. (6) of Ref. 24, which was derived by 
means of Eq. ( |l3[ ) , is a non-Hermitian matrix, and thus 
it is unreasonable to obtain real energy eigenvalues from 
this matrix. 

Figure 1 shows the E — <f> diagram obtained by diag- 
onalizing Eq. ( |l2] ) for the parameters t a = % = t c = 1, 
k = 0, and <f> — p/199 with 1 < p < 198 in the absence 
of the NNN hopping integral (t' a = t' b = t' c = 0). It is 
obvious that the energy spectrum is not symmetric about 
0=1/2 and there is no distinct array of energy eigenval- 
ues parallel to the energy axis in the two largest energy 
gaps, contrary to the argument of Gumbs and Fekete [see 
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FIG. 2. Energy eigenvalues versus <f> for a lattice with (a) 
(t a ,t b ,t c ) = (2,1,1) and (b) (t a ,t b ,t c ) = (4,1,1). Other pa- 
rameters are the same as Fig. 1. 

Fig. 1 and the relevant arguments of Ref. 24]. In fact, 
Fig. 1 of this paper is exactly the same as Fig. 3 of 
Ref. 21 presented by Claro and Wannier if E is replaced 
by — (E — eo)/2ei, which is quite natural because both 
of the two papers deal with the same problem. The only 
difference between the two papers lies in choosing the x 
and y directions; the axes taken in Ref. 21 are turned by 
7r/6 from those in this paper. Note that, even though the 
choice of the axes in Ref. 21 yields a three-term differ- 
ence equation [see Eq. (10) of Ref. 21] while the choice 
in this paper yields a five-term difference equation in the 
absence of NNN hopping integrals, the resultant energy 
band structures are the same, as shown above. 



A. Effects of anisotropic NN hopping integrals 



Figure 2 is the E — <j> diagram for the parameters 
tb = t c = 1, k = 0, and (a) t a — 2 and (b) t a = 4 in 
the absence of NNN hopping integrals. The increase of 
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total bandwidths and the occurrence of gap closing ex- 
cept for a few large ones are clearly seen. We argue that 
these behaviors may be generic effects of the hopping 
anisotropy, since the same-phenomena were also observed 
in the square latticeJla-LSNote that Fig. 3 of Ref. 24 is 
the E — <p diagram plotted for the same parameters used 
in plotting Fig. 2(a) of this paper, from which Gumbs 
and Fekete argued that the bottom of the energy band is 
very flat near E = —4 and that the spectrum is symmet- 
ric about = 1/2. However, it is evident that their argu- 
ments are incorrect. In regard to Figs. 2 — 4 of Ref. 24, 
we would like to point out that the E — (f> diagrams for 
the parameters (t a ,tb, t c ) — (1, 2, 1), (2, 1, 1), and (1, 1, 2) 
should be exactly the same when all the values of k in 
the magnetic Brillouin zone (MBZ; | fc^ | < 2n/Ma and 
\k y \ < 7r/\/3a) are taken into account, which means that 
there is no preferable direction in the lattice plane under 
a uniform perpendicular magnetic field. 
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FIG. 3. Rescaled energy eigenvalues as a function of t a for 
a lattice with (a) (j> = 1/5 and (b) (j> = 1/6, and tb = t c = 1- 
All the values of k in the MBZ are taken into account. 

In order to get more information on the anisotropy de- 
pendence of the subgap widths, we calculate the band 
structure as a function of the hopping anisotropy. Note 
that the energy eigenvalues in the absence of a magnetic 
field satisfy the condition \e(k)\ < 2(t a + t + t c ) [see 
Eq. (0)] and the total bandwidth increases linearly with 



the increase of t a - Assuming the situation also holds in 
the presence of a uniform magnetic field, we can intro- 
duce the rescaled energy defined by 



Er 



3E/{t a + t b + t c ). 



(15) 



Then, E re is always in the range of [—6, 6] regardless of 
the strength of the hopping anisotropy, and we can exam- 
ine more clearly the dependence of the subgap widths on 
the hopping anisotropy, in comparison with the isotropic 
case. 

Figure 3(a) shows the t a dependence of the band struc- 
ture for <j) — 1/5, where % = t c = 1 and all the values 
of k in the MBZ are taken into account. At isotropic 
case (t a — 1), there are five subbands and four subgaps. 
However, introducing the hopping anisotropy (t a =/= 1) 
changes the band structure considerably. Let us number 
the subgaps from the bottom of the energy spectrum. 
Then, in the regime of t a > 1, the first and the second 
subgaps close for a small amount of anisotropy, while the 
third and the fourth subgaps survive even for large val- 
ues of anisotropy. However, since the widths of the third 
and the fourth subgaps decrease with increasing t ai we 
expect that these subgaps will also close in the limit of 
t a — > oo. And, in the regime of t a < 1, the widths of 
the third and the fourth subgaps decrease with decreas- 
ing t a , as in the case of t a > 1- Meanwhile, the first 
and the second subgaps exhibit an interesting behavior. 
As t a decreases, the width of the first (second) subgap 
decreases such that it closes at t a ~ 0.8 (t a ~ 0.45). 
And, as t a decreases further, the subgap reopens and its 
width increases such that the width of the first (second) 
subgap becomes the same as that of the fourth (third) 
subgap, resulting in a symmetric band structure about 
E = in the limit of t a — > 0. Since the triangular lattice 
with the parameters (t a ,tb,t c ) — (0, 1, 1) is topologically 
equivalent to the square lattice with isotropic NN hop- 
ping integrals, it is natural of the band structure to be 
symmetric about E = 0. 

As an another example, we plot in Fig. 3(b) the t a de- 
pendence of the band structure for <f> — 1/6. At isotropic 
case (t a — 1), there might be six subbands and five sub- 
gaps. However, due to the occurrence of band touching 
between the second and the third subbands, the second 
subgap is missing and the energy spectrum consists of 
four subgaps. And, in the regime of t a > 1, we observe 
that the first subgap close even for a small amount of 
anisotropy and the third subgap close for a large value of 
anisotropy. Besides, since the widths of the fourth and 
the fifth subgaps decrease with increasing t a , the clos- 
ing of these subgaps is expected in the limit of t a — * oo. 
Meanwhile, the t a dependence of the missing second sub- 
gap exhibits an interesting feature; the subgap opens 
with increasing t a , even though its width is very small. 
In the regime of t a < 1 , we can see that the width of the 
first (second) subgap increases with decreasing t a such 
that the width becomes the same as that of the fifth 
(fourth) subgap, resulting in a symmetric band structure 
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about E = in the limit of t a — > 0. And, the width of 
the third subgap decreases with decreasing t a and closes 
at t a — 0. Thus, there occurs a degeneracy at the band 



center of the energy spectrum, as in the case of the square 
lattice with isotropic NN hopping integrals. 
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FIG. 4. Energy eigenvalues versus cf> for a lattice with (a) t' = 0.2, (b) t 1 = 0.5, (c) t' — 1.0, and (d) t' = 1.2, respectively. 
Other parameters are the same as Fig. 1. 



B. Effects of isotropic NNN hopping integrals 

In calculating energy eigenvalues with varying the 
NNN hopping integrals, we set t a = % = t c = 1 and 
assume that the NNN hopping integrals are isotropic, 
i.e., t' a = t' b = t' c (= t'). Figure 4 shows the E — <j) dia- 
gram for (a) t' = 0.2, (b) t' = 0.5, (c) = 1.0, and (d) 
t' = 1.2, respectively. We can see that the fractal struc- 
ture in the energy spectrum appears despite of introduc- 
ing the NNN hopping integral, which may attribute to 
the isotropy of the NN and the NNN hopping integrals. 
We can also see that, even though the total bandwidths 
for = and 1 increase monotonically with increasing t', 
the bandwidths for generic values of <f> depend crucially 
on t' , resulting in diverse shapes of E — <p diagrams. For 
small values of t' , the bandwidths become narrower than 
the case of t' = and the resultant diagram exhibits 
a vary narrow structure [see Fig. 4(a)]. However, the 
bandwidths increase with increasing t', resulting in a fat 
structure of E — 4> diagram [see Figs. 4(b) - 4(d)]. 

Figure 5 shows the t' dependence of the band struc- 



ture for <p = 1/2, which exhibits two critical values of 
t'. The subgap width increases with increasing t' up to 
t[ ~ 0.25, and decreases up to t' 2 = 1 such that gap 
closing occurs at t' 2 . And, then it increases again for 
t' > t' 2 . It may be interesting to compare this behavior 
with the case of the square lattice, since the triangu- 
lar lattice with t' ^ is topologically equivalent to the 
square lattice with isotropic (anisotropic) NNN hopping 
integrals if t' — 1 (t' ^ 1). In the square lattice, introduc- 
ing the NNN hopping integrals removes the degeneracy 
appearing at E — for even values of qS3 However, in 
the triangular lattice, a degeneracy occurs at E = in 
the energy spectrum for 0=1/2 when t' = t. 

The t' dependence of the band structures for = 1/3 
and 2/3 also exhibit interesting features. In these cases, 
there are three subgaps at t' = 0. And, with increasing 
t' , the widths of these subbands decrease up to a certain 
value of t' [see Fig. 4(a)], and then increase such that the 
two subgaps close [see Figs. 4(b)-4(d)] for large values of 
t'. Particularly, the lowest (highest) subband edge for 
0=1/3 (2/3) goes downward (upward) with increasing 
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FIG. 5. Energy eigenvalues as a function of t' for a lattice 
with = 1/2. All the values of k in the MBZ are taken into 
account. 
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FIG. 6. Energy eigenvalues versus for a lattice with t = 
and t — 1. The data are obtained for = p/199 with 
1 < V < 198 and for jfe = 0. 



t' such that £ ; (0 = 1/3) = E t ( 
E h ((j> = 0)] in the limit of t'/t - 



» = 1) [E h (<f> = 2/3) = 
oo, where Eu^) means 
the lowest (highest) subband edge. From this behavior, 
it can be expected that the periodicity of the E — di- 
agram in the limit of t'/t — > oo reduces by three times 
in 0. Figure 6 shows the E — <fi diagrams for the lattice 
with t = and t' = 1, which confirms this expectation. 
Note that it is nothing but an energy spectrum for the 
lattice with the isotropic hopping integral (t' = 1) and 
the lattice constant Via. In order to see the t' depen- 
dence of the band structure for generic values of <fr, we 
plot in Fig. 7 the E — t' diagrams for (a) <fi = 1/7 and 
(b) 4> — 2/7. We can see that gap closing, gap reopening, 
and band crossing occur depending on the values of t' . 

Before concluding this section, we discuss the symme- 
try of the energy spectrum. Denoting the set of energy 
eigenvalues for a wave vector k under a flux <f> as E(<j>; k) 
and the energy spectrum taking into account all the val- 
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FIG. 7. Energy eigenvalues as a function of t' for a lattice 
with (a) 0=1/7 and (b) = 2/7. All the values of k in the 
MBZ are taken into account. 

ues of k in the MBZ as E(<fi), we can see the follow- 
ing properties of the band structure for generic values 
of NN hopping integrals: (i) E(cj>;k) — E(cj>;—k), (ii) 
E(<f>) = E(-(f>), (iii) E{4>) = -E(l-(f>) with < < 1/2, 
and (iv) E{<j>) = £'(0 + 2). Note that the property (i) can 
be easily checked in Eq. ( |l2| ) ; A — * under the transfor- 
mation k — > —k. Note also that the property (ii) means 
that there is no way for an electron to discern the direc- 
tion of the magnetic field. Figures 1 and 2 show that the 
properties (iii) and (iv) hold in the presence of the NN 
hopping anisotropy. And, Fig. 4 also shows that the prop- 
erties (i)-(iv) hold in the presence of NNN hopping inte- 
grals except for the case of t = 0; when t = 0, the prop- 
erties (iii) and (iv) are replaced by E((f>) — —E(l/3 — 0) 
with < < 1/6 and E{<f>) = E((f> + 2/3), respectively, 
as can be seen in Fig. 6. 



IV. SUMMARY 

In summary, the effects of anisotropic NN hopping in- 
tegrals and of isotropic NNN hopping integrals on the 
magnetic subband structure of an electron on a triangu- 
lar lattice were studied within the tight-binding approxi- 
mation. It was found that the phenomena of bandwidth 
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broadening, gap closing, and gap opening are generic fea- 
ture of introducing the NN hopping anisotropy. It was 
also found that the magnetic subband structure changes 
considerably by introducing the NNN hopping integrals; 
bandwidth broadening, gap closing, gap reopening, and 
band crossing occur depending on the strength of the 
isotropic NNN hopping integrals. In the study, the de- 
pendence of the subgap widths on the hopping anisotropy 
and the NNN hopping integrals was illustrated in detail, 
and the symmetries of the magnetic subband structures 
with and without the hopping anisotropy and the NNN 
hopping integrals were also discussed. In order to deepen 
the understanding on the band structure of the triangu- 
lar lattice, a further study ^af the effect of, for example, 
nonuniform magnetic fieldsE3 and an interaction between 
electronscH on the magnetic subband structure of the lat- 
tice is required, and the work is on doing. 
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